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Abstract 

We present a generic model for the description of non-equilibrium Bose-Einstein condensates, 
suited for the modelling of non-resonantly pumped polariton condensates in a semiconductor 
microcavity. The excitation spectrum and scattering of the non-equilibrium condensate with a 
defect are discussed. 
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1. Model and elementary excitation spectrum 

Due to the finite polariton life time, the nature of polariton condensates [1] is fun- 
damentally different from their superfluid 4 He and atomic Bose-Einstein condensate 
counterparts: the polariton condensate has to be constantly replenished and arises as 
a dynamical equilibrium between pumping and decay. The theoretical modelling of such 
systems in which interactions, coherence, pumping and decay are equally important poses 
a challenge that was taken up only recently [2,3,4]. We will present here our mean field 
model of such condensates [4] and use it to study the excitation spectrum of a homo- 
geneous non-equilibrium condensate and the problem of a small defect moving through 
it. 

Our model does not include any details on the specific relaxation mechanisms of high- 
energy polaritons into the condensate, but is only based on some general assumptions: i) 
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a single state of lower polaritons is macroscopically occupied so that it can be described 
by a classical field; ii) The momentum space can be dcvided in two parts: one part 
at small momenta where the coherence is important and one part at large momenta, 
where it is negligble. The polaritons in the high-momentum states act as a reservoir that 
replenishes the condensate. Under typical excitation conditions, the wave vector scale to 
separate both systems should be chosen of the order of a few /im -1 ; iii) The state of the 
reservoir is fully determined by its spatial polariton density n R (x). This last asumption 
requires that the reservoir polariton momentum distribution reaches some stationary 
state in momentum space. Under these assumptions, the condensate dynamics is to a 
first approximation described by a generalized Gross-Pitaevskii equation including loss 
and amplification terms 
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[R(n R )- 1 ]+g\^\ 2 + 2gn R ^ 7 (1) 



where m is the lower polariton mass, 7 its decay rate and g is the strength of the polariton- 
polariton interaction within the condensate. The stimulated scattering of reservoir po- 
laritons into the condensate is modeled by the term R(n R ) and the mean field interaction 
experienced by the condensate polaritons due to elastic collisions with the reservoir po- 
laritons is given by 2gn R . The description (1) of the polariton condensate in terms of a 
deterministic classical field requires its density and phase fluctuations to be small. This 
regime is reached for pump powers well above the condensation threshold. The equation 
for the condensate dynamics is coupled to a diffusion equation for the reservoir polaritons 

= P- lR n R -R(n R )\^(x)\ 2 + DV 2 n R , (2) 

where P is the pump rate due to external laser, 7^ is the reservoir damping rate and D 
its diffusion constant. 

The stationary state and the elementary excitation spectrum of Eqns. (1) and (2) 
are discussed in Ref. [4]. In case the reservoir damping rate 7^ is much larger than the 
condensate polariton damping rate 7, the condensate excitation spectrum is of the form 



iT / r 2 

u±(k) = -— ± y uj B og{k) 2 - — , (3) 

Here, loboq = [(k 2 /2m + 2/z)(fc 2 /2m)] 1 / 2 is the usual Bogoliubov dispersion of dilute 
Bose gases at equilibrium. The non-equilibrium nature of the system is quantified by 
the effective relaxation rate T = £7, where £ depends on the pumping rate and on the 
functional form of R(n R ) [4]. The most important differences with the excitation spectrum 
of equilibrium condensates is the non vanishing imaginary part of oj(k) for all k ^ and 
the flatness of its real part for small k. The '+'-branch is diffusive for small wave vectors. 
A similar excitation spectrum was found in Ref. [2] for a specific model of nonequilibrium 
condensation, within a completely different approach, indicating that the form (3) is a 
general result for nonequilibrium condensates. 

Our model is also straightforwardly applied to the Joscphson oscillations between two 
condensates connected by quantum mechanical tunneling. The frequency of the density 
oscillations between the two wells is given by Eq.(3) with uusog replaced by the equilibrium 
Josephson frequency [8] . 
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2. Flow past a defect 



One of the benchmark properties of condensed Bose systems is superfluidity. As a 
first step in the study of superfluidity in non-equilibrium systems, we will discuss the 
scattering of a moving condensate on a defect [5]. Experimentally the condensate could 
be accelerated by applying an external force to it, e.g. by making a sample with a steep 
wedge in the cavity thickness or using surface acoustic waves to accelerate the polaritons 
[7]. Defects are naturally present in the form of disorder, but can also be deliberately 
created by structuring the cavity mirrors [6]. 

The perturbation on top of a condensate moving with velocity v due to a small defect 
potential Vdef — gS(r) at rest can be studied in perturbation theory [5]. The change in 
condensate wave function is in momentum space given by 

[^(k),^*(2k -k),5n r ] T = -C^ky^gMO), -<^o(«), 0] T , (4) 

where C v (k) the matrix from linearized motion equations (1),(2) around a steady state 
rii{(x) = no, ^(x) = °' x , where ko = mv. Much about the response of the flowing 
condensate is learned by studying the poles k res of L v (k) _1 in the complex k-plane, i.e. 

uj v ±(k res ) = uj±(k r es - k ) - v • (k res - k ) = 0, (5) 

We restrict now our attention to the perturbation of the wave function in the direction of 
the defect velocity. The main contribution comes from the Fourier component with k || v. 
The modulus of the resonant spatial frequency fc|| res is plotted in Fig.l. The panels show 
from left to right an equilibrium condensate (r = 0), a slightly non-equilibrium conden- 
sate (r//i = 0.01) and a condensate where interaction effects and losses are comparable 
(r = /i/2). Unlike for the temporal frequencies from Eq. (3), the imaginary part of the 
spatial frequencies should not be negative, but the spatial frequencies with a negative 
(positive) imaginary part correspond to positions on the left (right) hand side of the 
defect. 
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Fig. 1. The real (upper panels) and corresponding imaginary (lower panels) part of the spatial frequencies 
(in units of the inverse of the healing length £ = yj h/m^i) that are excited by a stationary defect in a 
condensate moving at a speed v for several values of T/fi. 

Let us start discussing the resonant spatial frequencies of the equilibrium condensate. 
The phase fluctuations, that give a pole at zero spatial frequency, do not couple to the 
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external potential. The two other branches of k\\ rea are purely imaginary for velocities 
smaller than the speed of sound, implying that the condensate perturbation is spatially 
damped. The excitations are only virtually excited by the defect and no condensate 
momentum is dissipated. At v = c, there is a bifurcation, where a pair of conjugate 
imaginary roots turn real, the excitations go 'on shell', are radiated and dissipate the 
condensate momentum. 

In the case of non-equilibrium condensates (central and right panels of Fig. 1), the 
bifurcation point is shifted to velocities v < c and Im(fcii res ) < 0. No spatial frequencies 
with zero imaginary part appear anymore. Therefore, in contrast to the equilibrium 
condensates, no sharp transition in the condensate perturbation as a function of its 
velocicy is expected. 



Fig. 2. The perturbation of the condensate moving condensate due to a small defect at x = 0. Condensate 
velocities: v/c = 0.7 (left panel), v/c = 1.3 (center panel) and v/c = 2 (right panel). Effective damping 
rate T = /j,/2. 

This is confirmed by the condensate perturbation in real space, shown in Fig. 2 for 
several velocities v. The condensate perturbation depends sensitively on the velocity v, 
but no sharp transition occurs. It is smoothened due to the finite polariton life time. 
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